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We consider one-loop self-energy and vacuum polarization radiative corrections to the shift of 
atomic energy level due to finite nuclear size. Analytic expressions for vacuum polarization correc- 
tions are derived. For the self-energy of p\/2 and P3/2 states in addition to already known terms 
we derive next-to-leading nonlogarithmic Za-terms. Together with contributions obtained earlier 
the terms derived in the present work give explicit analytic expressions for Si/j and P\j2 corrections 
which agree with results of previous numerical calculations up to Z=100 (Z is the nuclear charge 
number). We also show that the finite nuclear size radiative correction for a pz/2 state is not small 
compared to the similar correction for a P\/2 state at least for small Z. 
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I. INTRODUCTION 

Experimental and theoretical investigations of the radiative shift (Lamb shift) of energy levels in heavy atoms and 
ions is an important way to test Quantum Electrodynamics in the presence of a strong external electric field. One 
of the effects related to this problem is the dependence of the Lamb shift on the finite nuclear size (FNS). One can 
also look at this effect from another point of view. It is well known that there is an isotope shift of atomic levels due 
to the FNS. Corrections we are talking about are the radiative corrections to the isotope shift. There are two types 
of corrections. The first type is due to vacuum polarization and we will call them VPFNS corrections. The second 
type is due to the electron self-energy and vertex and we will call them SEVFNS corrections. Calculations of both 
VPFNS and SEVFNS corrections have attracted a lot of attention, for a review see Ref. 0. Until recently it was 
mainly numerical work. The corrections for Is, 2s, and 2p states have been calculated numerically, exactly in Za, 
in Refs.0, 00,0, 0,013; see also review 0] for muonic atoms . Here Z is the nuclear charge number and a is the 
fine structure constant. Analytical calculations have been based on Zct-expansion. The s-wave VPFNS correction 
has been calculated a long time ago in order a(Za) 0,0]. An ultraviolet logarithmic enhancement of higher order 
contributions to the VPFNS correction has been revealed in Ref. 0]. The enhancement factor is \n 2 (Xc/ro), where 
Ac is the Compton wave length and ro is the nuclear radius. In Ref. |12| the double-logarithmic contribution to the 
radiative correction to atomic parity nonconservation has been calculated analytically exactly in Za. The formula 
derived in 0] gives also one of the double- logarithmic contributions to the VPFNS correction for Si/ 2 and Pi /? states. 
The SEVFNS correction for an s-wave state has been calculated analytically in order a(Za) in 

Refs.millllllH, 

see also reviews ITtI fl8l. Structure of higher-order contributions to the SEVFNS correction has been elucidated in 
our recent papers |1 all fJ|. The contributions are also logarithmically enhanced, but the enhancement factor contains 
only first power of ln(Ac/?"o)- For Si/ 2 and pi/ 2 states prefactor before the logarithm has been calculated in Ref. 0] 
in order a(Za) 2 and in Ref. 0] exactly in Za. Finally, contributions of the order a(Za)° to SEVFNS corrections 
for P1/2 and p 3 / 2 states have been calculated in recent papers p^.l2fj|. 

In the present work we first calculate nonlogarithmic a(Za) terms of SEVFNS corrections for p\j 2 and P3/2 states. 
Together with previously known contributions this gives a full analytic description for P1/2 SEVFNS correction which 
agrees with previous numerical results up to Z = 100. Our calculation demonstrates also that the p 3 / 2 SEVFNS 
correction is comparable with the P\j 2 SEVFNS correction. The second part of the work is devoted to analytical 
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calculations of VPFNS corrections for s\/2 and p\/2 states. Analytic formulas obtained in this paper agree with 
previous numerical results up to Z = 100. 

The structure of the present paper is as follows. In Sec. II the general structure of SEVFNS and VPFNS corrections 
is discussed. In Sec. Ill we calculate a(Za) terms for P1/2 and p 3 / 2 SEVFNS corrections, summarize all known terms 
for s-wave and p-wave corrections and compare analytic expressions with previous numerical results. In Sec. IV we 
first calculate unknown contributions to the VPFNS correction. Then we summarize contributions calculated in the 
present work and corrections obtained before and compare our analytic results with previous numerical calculations. 
Finally, Sec. V presents our conclusions. 



II. GENERAL STRUCTURE OF SEVFNS AND VPFNS RADIATIVE CORRECTIONS 



Due to the finite nuclear size, the electric potential V(r) of the nucleus differs from that for a point-like nucleus. 
The deviation is 



SV(r) = V{r) - 




Throughout the paper we set H = c = 1. The diagram that describe the FNS effect in the leading order is shown 
in Fig^,. The double line corresponds to the exact electron wave function in the Coulomb field, the zigzag line 





d e 

FIG. 1: The leading contribution to the FNS effect is given by diagram (a), and one loop radiative corrections to the effect are 
given by diagrams (b-e). The double line denotes exact electron Green's function in the Coulomb field of the nucleus, the cross 
denotes the nucleus, the zigzag line denotes the FNS perturbation [7]), and the dashed line denotes the photon. 



with cross denotes the perturbation l|T|l. and the dashed line denotes the photon. In momentum representation the 
perturbation reads 

6V(k) = ^[l-F(k*)], (2) 

where ^F{k 2 ) is the nuclear form factor. An important fact is that is a local interaction which is nonzero only at 
r ~ tq, where tq is the nuclear radius. For estimates one can use the following formula for tq 



1.1 A 1/3 fm w 1.5 Z 1/3 fm , (3) 



where A is the mass number of the nucleus. Dynamics of electrons at r ~ ro is ultrarelativistic at any value of Z. 
The Dirac wave function at r <C ZaXc is of the form, see Ref. 21] 

m = Nr^^-^) , (4) 
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where SI and ft = — (a ■ n)il are spherical spinors; k = —1 for s\/ 2 state, n — 1 for p 1 / 2 state, and n — —2 for p 3 / 2 
state; 7 = n 2 — (Za) 2 ; and N is a constant known for each particular state. For hydrogenlike ions one can find 
N using wave functions presented in Ref [2l]]. For multi-electron ions or atoms the constant N can be calculated 
numerically. It can also be calculated analytically but in the semiclassical approximation 22] . 

Due to ultrarelativistic nature of the problem the strong relativistic enhancement is a special property of SEVFNS 
and VPFNS radiative corrections. This is similar to atomic parity nonconservation (PNC) [22j . The relativistic 
enhancement factor for Sx/2 and V±l2 states is proportional to R re i ~ (Ac /Zaro) 2 ^ 1 ^^, The factor is R re i ~3 for Cs 
and R re i « 9 for Tl, Pb, and Bi 22]. The enhancement factor R re i is divergent at tq — > 0. As a result behavior of 
SEVFNS, VPFNS, and PNC radiative corrections differs essentially from that of radiative corrections to the hyperfine 
structure. 

The FNS energy shift AE in the leading approximation, diagram Fig. la, has been calculated in numerous works 
both analytically and numerically. For a hydrogenlike ion there are simple parameterizations suggested in Ref. [23j 



Z 2 , , , 2 J2ZuR\ 21 



Z* a 2 (n 2 -1) M (2ZaR\ 2 ^ 

The energies are given in atomic units, E = ma 2 = 27.2eV; m is the electron mass; n is the principal quantum 
number, n = 1, 2, ...; R ss r$ is some effective radius (for details see Ref. |23|); and coefficients fi — fi(Za) are 

f 1Sl/2 (Za) = 1.380 - 0.162(Za) + 1.612(Za) 2 , 

f 2si/2 (Za) = 1.508 + 0.215(Za) + 1.332(Za) 2 , 

f 2pi/2 {Za) = 1.615 + 4.319(Za)-9.152(Za) 2 + 11. 87(Za) 3 . (6) 

To be absolutely correct we should say that expectation values of the perturbation Q taken with the wave functions 
are not sufficient to calculate the energy shifts. One has also to take into account deformation of the wave function 
inside the nucleus. However, Eqs. (0) take this effect into account, and as soon as it is done there is no need to care 
about this detail any more. 

Diagrams Fig^b-e) show the radiative corrections under discussion, see also comment |24| . Diagrams FigOJ) (self- 
energy) and Fig^t (vertex) describe the SEVFNS correction and diagrams Fig^i and Fig^3 describe the VPFNS 
correction. Note that Fig^i corresponds to a modification of SV (see Eq. (Q) due to the vacuum polarization, 
and Fig^ corresponds to a modification of the electron wave function due to the vacuum polarization. Each of the 
diagrams FigQJb-e) gives some energy shift 5Ei, i=b,c,d,e. Let us stress that for S\i 2 and P\i 2 states we always 
calculate a relative value of the radiative correction 

A Sl/2 = SE Sl/2 /AE Sl/2 , 

A Pl/2 = SE Pl/2 /AE Pl/2 , (7) 

where AE is the FNS shift given by diagram FigQJi for the same atomic state. To find the value of AE one can 
use Eqs. |J5J, or a numerical calculation or any other method. Situation with p 3 / 2 correction is different and we 
discuss this at the end of section III (see Eq. Ij23(l . So, a crucially important point is that we calculate analyti cally 
only the relative quantity A . This is the central issue of this work as well as of our recent papers 0, 0, llii Il9| . 
The point is that convergence of Za-expansion for A is quite reasonable while the ultrarelativistic character of the 
FNS effect mentioned above makes the convergence of Za-expansion for 5E very poor. The good convergence of the 
series for the quantity A is related to separation of scales. Radiative corrections which do not contain ln(Ac/'*o) are 
related to quantum fluctuations at "large" distances, r ~ Ac- On the other hand, the perturbation is nonzero only 
at small distances r ~ tq <C Ac- Therefore, the enhancement related to small distances in the contribution of "large 
distance" quantum fluctuations is factorized similar to (J5J , in essence this is a kind of effective operator approach. 
There are also quantum fluctuations from distances r < r < Xc, their contribution cannot be trivially factorized 
and we take special care about these fluctuations. This is where terms dependent on ln(Ac/Vo) come from. The 
factorization of small distance enhancement in the part of the FNS effect related to the vacuum polarization was 
implicitly demonstrated in Ref. 25]. 



III. SELF-ENERGY AND VERTEX FNS RADIATIVE CORRECTIONS 



Technically the most complicated are the self-energy and the vertex FNS radiative corrections (SEVFNS) given by 
diagrams in Fig^) and Fig[^. It has been demonstrated in paper ^^] that relative SEVFNS corrections for si/ 2 and 
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Pi/2 states are of the form 



AHb\ c /r )+B. 



(8) 



where A and B are functions of Za independent of ro; b = exp(l/(27) — C — 5/6), 7 = wl — (Za) 2 , C w 0.557 is the 
Euler constant. The function A(Za) is the same for Sx/2 and P1/2 states. This function has been calculated in the 
leading order in Za in Ref. [15j, the result reads 



A = - 



a{ZaY 



15 7T Z 



(9) 



It has been also calculated exactly in Za in Ref. ^(| • Plot of A from Ref. 0] is presented in FigEl Deviation of © 



s 

N 




100 



FIG. 2: The function A , Eq. calculated to all orders in Za. Value of the function is given in units a(Za) 2 /tv. 



from exact A(Za) is small even at Z = 100. The function B(Za) for an S1/2 state has been calculated in the leading 
order in Za in Refs. [la, LLal 
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Z? s = -a(2a) [ — -4 In 2 



(10) 



It agrees numerically with earlier results 0,0]. Combining Eqs. (|5)l.(|9")l. (|10[l one gets the the following expression 
for s-wave SEVFNS correction 



aSEV 

Sl/2 



v ^ ,'23 „, \ a(Zaf ( 15 ir 2 \ 



ln(6Ac/ro) 



(11) 



The leading unaccounted term is of the order ~ a(Za) 2 /ir. The correction (|ll|l in % is plotted in Fig|31by the dashed 
line. The same correction, but with A taken from Figgis plotted in Fig|3|by the dotted line. Our analytical results 
for Afj^ are in a good agreement with previous numerical data 0, |{| shown by squares and stars for Is state and 

by triangles and circles for 2s state. Some dependence of Af ^ on the principle quantum number n can appear from 
function B s in the order a(Za) 2 /n. However, one can see from the numerical data 0,0 presented in Fig. [3]that this 
dependence remains very weak up to Z = 100. 

The function B p i/ 2 for p 1 / 2 states in the leading order has been calculated in Refs. flU |20|. 



B 



(0) 
pl/2 



a — I In -T7 

9tt V (Za) 2 



(12) 



For 2p-| /2 state the constant is d = 0.910. This constant depends slightly on the principle quantum number n, see 
Ref. [Si. 

It is interesting to note that Zen-expansion of B v \i% starts from ahx(Za) term, while the same expansion for B s 
starts from a(Za). It is related to a different infrared behavior of quantum fluctuations, see discussion in Ref. 19] . 
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FIG. 3: Relative SEVFNS corrections (%) for S1/2 and pi/2 states. The dashed line shows the correction A s given by Eq. 
11 1\) . The same correction, but with A taken from Fig\^ is plotted by the dotted line. Results of computations of A?^ for 
Is states are shown by squares and stars 0/. Similar results for 2s states are shown by triangles Q/ and circles 0/. The 
solid line shows the correction Ap®^ given by Eqs. (51). [PP . 121$ . The same correction, but with A taken from Fig\Q is plotted 

by the dashed-dotted line. Results of computations of A^^ for 2p 1 / 2 states are shown by diamonds 0/ and crosses The 
correction Ap^ given by Eq. \24\ is plotted by the long-dashed line. 



Now we calculate the a(Za) term in the expansion of B p i/2- The electron wave function at small distances is given 
by Eq. |0J. Let us first calculate the leading contribution to the FNS shift given by diagram in Fig^,. For Si/ 2 
and p 3 /2 states the upper component of the Dirac spinor is much larger than the lower one. Hence, the upper 
component determines the FNS shift of such a state. On the other hand, for a pi/2 state the lower component and 
hence its contribution to the FNS shift is dominating. A straightforward calculation gives the following values for the 
FNS shifts of s and p states (diagram Fig^a)) 

AE nSl/2 = A( Za )4 m 3 < r 2 > ) 

n 2 — 1 

AE npi/2 = -^(Za) 6 m 3 <r 2 >, 
n 2 - 1 

AE np3/2 = —^(Zafm 5 < r 4 > . (13) 

Here < r 2 > and < r 4 > are values of r 2 and r 4 averaged over charge density of the nucleus, n is the principle 
quantum number. 

The low-momentum expansion of the nuclear electric form factor, see Eq. J2J, is of the form 

Hk 2 )*l-^ <r 2 >+^<r 4 > ■ (14) 

Modeling the nucleus as a uniformly charged ball one gets 

< r 2 >= I rg, < r 4 >= | r\ . (15) 

As one should expect, the FNS corrections (|T3|l obey the following inequalities AE Sl/2 ^> AE Pl/2 ^> AE P3/2 . 

We consider now radiative corrections which originate from quantum fluctuations at distances r > Ac- To calculate 
these "soft" corrections it is sufficient to use nonrelativistic electron wave functions and, instead of 0]), to use an 
effective FNS perturbation that reproduces the FNS correction AE nsl / 2 for s-wave states, see Eq. (f*H?|) and Ref. 19]. 
The effective perturbation reads 

5V eff (r)=gS(r) , .9 = ^ < r 2 > . (16) 
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In the leading approximatio n SE VFNS corrections for I'pxii and 2p 3 / 2 states, corresponding to diagrams Fig^bjC), 
have been calculated in Refs. 19, 20] : 



SE {0) - 
where F 



1 



(Za) 2 
1 



(Za) 2 
a(Za) 5 gm 3 
48^ 



0.910 
0.215 



(17) 



These corrections are due to quantum fluctuations with the frequency (Za) 2 m < ui < m. Therefore, the effective 
operator is calculated with the relativistic effects taken into account. However, due to the behavior of the wave 
functions of p-states at r <C Xc/(Za), it is possible to average this operator in the leading approximation over the 
nonrelativistic wave function. For s-states the situation is quite different |26| . 

The ratio 5E$ /AE 2pi/2 defined by Eqs. |T7|) and gives the relative correction (|12fl . This means that Eqs. 



(1 1 71) correspond to zero order in (^a)-expansion of the function B Pl 



We will see that next-order corrections SEp]^ /2 and SE P ^' do not contain any logarithms, hence they come from 
quantum fluctuations at r ~ Ac- Therefore it is convenient to calculate these corrections using the effective operator 
approach, see, e. g. Ref. 01- ^ n this approach we first of all calculate scattering amplitude, described by diagrams 
shown in Fig0] External electron momenta in Fig0] are on mass shell and it is sufficient to consider scattering at 
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FIG. 4: Linear in Za contributions to the electron scattering amplitude related to the p-wave SEVFNS correction. The solid 
line describes electron, the zigzag line denotes FNS perturbation 11 bV . the dashed line denotes photon, and the cross denotes the 
nucleus. The diagrams with the permutation of Coulomb and perturbation lines should be added. 

low momenta \pi^\ <C m. Here p\ and p 2 are initial and final momenta of the electron, respectively. To find the 
energy shift one has to average the scattering amplitude over nonrelataivistic wave function of atomic electron. We 
have already used this technique to calculate l|10|) and to calculate a similar term for atomic parity nonconservation 
effect, see Refs. Similar to these papers we use here the Fried- Yennie gauge [27]]. One can find expressions 

for renormalized self-energy and vertex operators (ultraviolet renormalization) in this gauge in Refs. ^3] and |l6l. 
However, the present calculation is more complicated than that in [Til Utij . The point is that in Refs. 0, lid] it 
was sufficient to calculate only the forward scattering amplitude, p\ = p 2 . In the present case we have to consider 
scattering at arbitrary angle, p\ ^ p 2 . It is not just a simple technical complication. The forward scattering 
amplitude is always infrared convergent. On the other hand the finite angle amplitude is always infrared divergent 
because of long range nature of the Coulomb field. To solve the problem we regularize the Coulomb interaction 
in Fig0J 1/fe 2 — » l/(k 2 + A 2 ), where m 3> A ^> |pi,2|- Performing calculations we throw away all terms inversely 
proportional to A because they correspond not to the true radiative corrections, but just to Coulomb corrections to 
the scattering amplitude corresponding to Eqs. I|17|) . Since we perform calculations in the Fried- Yennie gauge there 
is no need in an infrared regularization of non-Coulomb photon propagator in diagrams FigU After a pretty long 
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calculation we have found the following expressions for diagrams in Fig^ 



M a = G 



Mh + M c = G 



M, = G 



— (cr ■ p 2 )(cr ■ px) 



1 



s 



227 

(cr • p 2 )(cr ■ px) + ( — 



Gl 



--+ln2 ) (cr ■ p 2 )(cr ■ Pl ) - w+ o ln2 (Pi-Pa) 



iln2 



(Pl -P2^ 



288 



(18) 



Here <x is the Pauli matrix, and G = ga(Za)/m 2 . Thus, the total scattering amplitude corresponding to Fig^reads 



M = M a + Mb + M c + Md = G 







-(Jr 









559 

(cr -p 2 )(cr -px) + ( — -31n2 ) (px -p 2 ) 



(19) 



To find the effective Hamiltonian we transfer (|19l) to coordinate space using substitutions (cr ■ p 2 )(cr ■ px) — > (tr • 
p)S(r)(cr ■ p), and (pi • p 2 ) — + (p<5(f)p) , where p = — iV. Finally, calculating expectation values of the effective 
Hamiltonian with the wave functions of np states we obtain the following corrections to energy levels : 



/l ^ -y ln2 

r L 4 In 2 

1 432 



a(Za)V 3 (n 2 -l) 
Ann 5 



Using tfi^ 2 from Eq. lfP7jl. 5E$ 1/2 from Eq. J20J), AE npi/2 from (O as well as definitions Q, 



expression for B p x/ 2 valid up to the first order in (Za). 

1 



15 



pl/2 



97 1 ln 



(Zaf 



0.910 



, 379 16, 



(20) 

we find explicit 
(21) 



Strictly speaking the constant near the logarithm depends on the principal quantum number. Equation l|21|) corre- 
sponds to the 2px/2 state, for the 3px/ 2 state one has to replace 0.910 — > 0.908, see Ref. [2JJ. Combining Eqs. JHJ,© 
and (|21|l one gets the the following expression for the SEVFNS radiative correction for a px/ 2 state 



A 



SEV 
Pl/2 



9tt 



In- 



(Zaf 



0.910 



a{Za) 



379 
432 



1G 



ln2 



a{Za) 



15 

T 



7T" 

"g 



ln(6A c /r ) . (22) 



The leading unaccounted term is of the order ~ a(Za) 2 /%. The correction 11221) in % is plotted in FigOby the solid 
line. The same correction, but with A taken from Figgis plotted in Fig|3]by the dashed-dotted line. Our analytical 
calculation of A p ^ agrees very well with previous numerical results for 2px/ 2 states shown by diamonds jj] and 
crosses see also comment p8|. 

Let us discuss now the SEVFNS correction for a p 3 / 2 state. In this case the definition (7J for the relative value of 
the correction is not sensible. The point is that according to (|17f) and l|20|> the SEVFNS correction SE P3/2 oc (r 2 ), and 
according to l|13|) the FNS correction AE P3/ ^ oc (r 4 ). Therefore the definition J7J) would imply that A P3/2 ^> 1. This 
itself is not a problem, just some rescaling. The real problem is that different powers of r 2 indicate that relativistic 
factors for 5E P3/2 and AE P3/2 are different. Therefore convergence of (Za)-expansion for A P3/2 defined according to 
Eq. must be very poor. Instead we define A P3/2 as 



^ n 2?3/2 — ^Fj np3/2 / AE npi/2 



(23) 



This is the ratio of the SEVFNS energy shift for np 3 / 2 state and the FNS energy shift for npx/ 2 state. Using definition 
(|23|l we absorb all small-distance (r ~ ro) effects into AE npi/2 , and therefore, according to the effective operator 

approach explained in section II, we expect a reasonable convergence of (Za)-expansion for A np3/2 . 



from Eq. l(T7|) . SE„J 3/2 from Eq. (t2U|) . AE npi/2 from (|TSf> as well as the definition we obtain 



Using SEi°J , 



A SEV 

P3/2 



(Za)' 



- 0.215 



7 I 559 

a(Za) 

V M 432 



4 In 2 



(24) 
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Again, the constant near the logarithm depends on the principal quantum number. Equation l|24[) corresponds to the 
2^3/2 state, for the 3p 3 / 2 state one has to replace —0.215 — * —0.217, see Ref. [2(|. There are no ln(Ac/r ) terms in 
A P3/2 , so the leading unaccounted term in (|24|l is of the order ~ a(Za) 2 /%. The correction A p ^ {%) given by Eq. 
(l2"4"j) is plotted in FigEUby the long-dashed line. Note that it is comparable to £Xjj,®Y ■ It means that absolute SEVFNS 
energy shifts SE P3/2 and SE Pl/2 are also comparable. 



IV. VACUUM POLARIZATION FNS RADIATIVE CORRECTIONS 



In this section we consider vacuum polarization FNS radiative corrections (VPFNS) given by diagrams in FigCJi 
and Fig^. Let us start from sj/2 states. In the leading order in (Za) the diagrams FigQi an d Fig^ can be 
represented as it is shown in FigEl Calculation of these diagrams for forward scattering p\ = p 2 is very simple, they 




FIG. 5: Linear in Za FNS vacuum polarization corrections corresponding to Za-expansion of diagrams Fig^l and Fig^jp. 
The solid line describes electron, the zigzag line denotes FNS perturbation 116]) . the dashed line denotes photon, and the cross 
denotes the nucleus. 

give equal contributions, and their total contribution to the s-wave VPFNS correction is |0, 

= \a{Za) . (25) 
Logarithmically enhanced contribution of diagram Fig^ has been calculated in Ref. 01 j it reads 



3vr7 \ ro 

This formula is equally valid for S1/2 and pi/2 states, for definition of parameters see Eq.JHJ). A logarithmically 
enhanced contribution of the diagram Fig^i has not been calculated yet. A similar contribution has been estimated 
in Refs.pjlll^| in relation to atomic parity nonconservation (PNC). For PNC this contribution is strongly suppressed 
by the factor 1 — 4sin 2 #n/, where 9\y is the Weinberg angle, therefore there have been no need in an accurate 
calculation. Now we perform an accurate calculation, the calculation is equally valid for s\/2 and P1/2 states. 
According to 10} the electron density at r <C (Za)Xc is 

p{r) = Mr 2( ~>-^ , (27) 

where M is some constant. Fourier transform of the density reads 

p(r)dr = sin (tt 7 ) T(2 7 ) , (28) 

where T(x) is the gamma function. Then, due to Eqs. (Q) and © the FNS energy shift given by diagram Fig^ reads 

AE = [ SV{r)p(r)dr = MAZasin (7r 7 )r(2 7 )ro 7 / -^-[1 - F\{y)] • (29) 

J Jo y 1 

where y — k 2 r 2 , and T\(y) = T(yjr^). 

To generate diagram Fig^i one has to account for the electron loop in the interaction energy 129fl . To do so we 
perform the following replacement in (|29(l . see, e.g., Ref. ^3] and comment (2^| 

1 a f 1 v 2 (l-v 2 /3) 

dv TTT, /„2 ' ( 3 °) 



k 2 nm 2 J 4 + (1 — v 2 )k 2 /m 2 
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This leads to the following expression for diagram FigCJi 



SEW = MAZa sin (7r7)r(27)r 2 



2T 



dy 



o y 



1+ 



rti-^ifo)] - 

1 TT 



1 dvv 2 (l-v 2 /3)(y/a 2 ) 
4+(l-v 2 )(y/a 2 ) 



(31) 



where a — mro = tq/Xc <C 1. Having in mind this inequality the expression i|31|l can be transformed to 
6E& = M4Zasm(n 7 )r(2 7 )rp(^)y™ 



ln^-- 
a 2 3 



where ^(0) = ^ 



3^(0) q 2(1 ~ 7) 7 r 3 / 2 4 1 -^r( 7 + i) 

+ 4( 7 - 1) sin(7r 7 ) r( 7 + 3/2) 
Finally, introducing notation 

dy 



J = 



(32) 



(33) 



and substituting Eqs. (|32|l and l|29|) in the definition {7J) we obtain the following expression for logarithmically 
enhanced contribution to the relative correction corresponding to Fig^i 



21n- 



5 d In J 

3 ~ <9 7 



3.F{(0)(AcA ) 2 ( 1 -^ tt 3 /^ 1 ^ r( 7 + 1) 
4( 7 - 1) sin(7T7) r(7 + 3/2) J 



(34) 



This formula is equally valid for s\/2 and pi/2 states. 

For large Z, as soon as (1 — 7) ln(Ac/?"ci) ^> 1, the .F-J-term in 113 111 can be neglected and 



A (d,log) _^ (<*_ 

V3tt 



A c 5 d In J 
2 In — 

r 3 d7 



(35) 



A similar high Z estimate has been used in |0, For small Z, more precisely for (1 — 7) ln(Ac/7"o) <^ 1, the 
correction behaves as 



a(Za) 
3tt 



(36) 



The simplest way to obtain this limiting behavior is to expand intermediate expressions (|29(l and (|31|) . however, one 
can also obtain the same result expanding the exact answer Ij34|) . 

Behavior of the exact expression (|34|l between asymptotical regimes (|35|) and <|36[l depends on the particular form 
factor F. However, the dependence is pretty weak, and for further analysis we use the form factor of uniformly charged 
ball of radius rp. 



sin^/y 



An explicit integration in (|33|l gives 



J 



Vv 

6cos(7r 7 )r(-l - 2 7 ) 



cos^y 



(37) 



d In J 



3 + 2 7 
7r tan(7r 7 ) 



+ 2^»(-l-27), 



(38) 



9 7 v " " ' 3 + 27 

where ip(x) = dh\Y(x)/dx. Combining together Eqs. i|25|) . i|26|l . (|34|l . I|38(l . as well as unaccounted in the present cal- 
culation correction u(Za) ~ a(Za) 2 /ir we obtain the following expression for the total vacuum polarization correction 



for s-clectron 



VP 

^1/2 



-a(Za) 



2a(Za) 2 



( a 



2 In 



37T7 

A c 5 
r Q 3 



In 2 



b\ c 



7rtan(7T7) 



3 + 2 7 



+ 2V>(-l-2 7 ) 



r 3 / 2 4( 1 -T)(3 + 2 7 )r( 7 + l) 



40sin(27T7)(7 - l)r(-l - 2 7 )r(7 + 3/2) V r 



c 



2(7-1)" 



■ u(Za) 



(39) 
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The correction <|39[1 with u(Za) = is plotted in Fig|HJi by the solid line. Results of computations of the same 
correction performed in Ref |3j] for Is states are shown by squares and for 2s states are shown by triangles. Agreement 
between our analytical calculation and Q) is pretty good. Moreover, we can determine u{Za) fitting the numerical 
data. Result of the fit reads 

a(Za) 2 , 
u = 2.5— . (40) 

The correction <|39|) with u(Za) given by 140|) is plotted in Fig|S^ by the dashed line. 





FIG. 6: Relative VPFNS corrections (%) for s-i/2 and pi/2 states. The solid line in Fig. "a" shows the correction " given 
by Eq. 1396 with u(Za) — 0. The dashed line shows the same correction with u(Za) given by Eq. Results of computations 

0/ are shown by squares for Is states and by triangles for 2s states. The solid line in Fig. "b" shows the correction given 
by Eq. [4 St with u(Za) — 0. The dashed line shows the same correction with u(Za) given by Eq. \4tfy - Results of computations 
for 2p!/ 2 states are shown by diamonds. 



Let us consider now the VPFNS correction for api/2 state. Logarithmically enhanced contributions (I26|l and l|34|) 
are exactly the same as for Si/ 2 states. There is also a contribution of zero order in (Za) calculated in |19| 



aVP(O) 

Pl/2 



45-7T 



(41) 



To calculate the first order correction we once more employ the effective operator approach and consider the scattering 
amplitude described by Fig|Sl in this case pi ^ P2- Similar to the procedure described in Section III we regularize the 
Coulomb interaction in Fig|SJ 1/fe 2 — > l/(fc 2 + A 2 ), where m A ^> |pi a|- Performing calculations we throw away 
all terms inversely proportional to A because they correspond not to true radiative corrections, but just to Coulomb 
corrections to scattering amplitude corresponding to Eq. I|41|) . Then we obtain the following result for the first-order 
correction 



pi/2 72 



(42) 



The coefficient is 23/72 = 1/8 + 7/36, where the first term comes from the diagram Fig[S^, and the second term comes 
from the diagram FigEJi. Combining together Eqs. (|41|l . (|42Jl . l|2f)JI . (|34|l . (|38Jl . as well as unaccounted in the present 
calculation correction u(Za) ~ a(Za) 2 /tt we obtain the following expression for total vacuum polarization correction 
for a pi/2 state 



a vp 

Pl/2 



-—a+-a(Za) 



(3 



21n^ 
3 



2a(Za) 2 nx c \ 
37T7 V r o / 

f 7r tan(7T7) 



r 6 ■ " 3 + 2 7 

^3/2 4 (i- 7 ) ( 3 + 2 7 )r( 7+ l) /Ac 
40sin(27r 7 )(7 - l)r(-l - 2 7 )r( 7 + 3/2) ^ r 



+ 2^(-l-2 7 ) 



2(7-1)' 



+ u(Za) . 



(43) 
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In principle the unaccounted contribution u(Za) can be different from that in l|39|) . The correction 1|43|) with u(Za) = 
is plotted in FigEJ) by a solid line. Results of computations of performed in Ref for 2p 1 / 2 states are shown 

by diamonds. Once more, the agreement is pretty good. We also find unknown u(Za) fitting the numerical data 0. 
The fit gives the result very close to that for s-wave, see Eq. (|40fl . The dashed line in FigJHJa shows the correction 
g3J with u(Za) given by Eq. I|40|). 

Finally we consider the VPFNS correction for &p 3 / 2 state^the relative correction is defined according to Eq. 1|23[) . 
There is a contribution of zero order in (Za) calculated in |l9| 

A^ p (°) = -—a . (44) 

It is equal to that for a pi/ 2 state, see Eq. I|41|l . The first-order correction is given by diagrams in FigEl The 
calculation of the correction, which is similar to the calculation of 1421) , gives 

= T2 a{Za) ■ (45) 

This correction comes from the diagram Fig|SJi only. There are no logarithmically enhanced contributions to Aj 3 P , 
therefore with accuracy ~ a(Za) 2 /ir accepted in the present work the correction A^ p is given by sum of 144(1 and 

^ = ~^ a + T2 a{Za) - (46) 



V. SUMMARY AND CONCLUSIONS 



In the present work we have calculated analytically the finite nuclear size effect on Lamb shift of S1/2, Pi/2; an d 
P3/2 atomic states. Corresponding radiative corrections are strongly affected by ultrarelativistic behavior of electrons 
in the vicinity of the nucleus. As a result of this behavior the convergence of Za-expansion is very poor for absolute 
energy shifts, and on the other hand the convergence is reasonably good for relative corrections which we consider. 
We define the relative corrections for s\/ 2 and pi/ 2 states according to Eq. (Q, and for p 3 / 2 states according to Eq. 
(I23|l . The self-energy and vertex finite size (SEVFNS) relative radiative correction for a S\/ 2 state is given by Eq. 
(fTTfl - The SEVFNS correction for a p 1 / 2 state is given by Eq. J22J), and the SEVFNS correction for a p 3 / 2 state is 
given by Eq. I)24[l. The vacuum polarization finite size (VPFNS) relative radiative correction for a S\/ 2 state is given 
by Eqs. ((231,(1301). The VPFNS correction for a p 1/2 state is given by Eqs. (gSJ , GUI) , and the VPFNS correction 
for a p 3 / 2 is given by Eq. (|46|) . Our analytical results for si/ 2 and p\/ 2 states agree well with previous numerical 
calculations. 

Finally, in Fig. 0we plot the total radiative corrections which include self-energy, vertex and vacuum polarization 
diagrams (SEVFNS + VPFNS). For s\/ 2 andpi/2 states at large Z there are strong compensations between SEVFNS 
and VPFNS contributions. We would like to attract attention to the p 3 / 2 correction. The same normalization is used 
for p\j 2 and p 3 / 2 states. Therefore, from Fig. \7\we come to a somewhat paradoxical conclusion that at large Z the 
finite nuclear size effect on Lamb shift of p 3 / 2 states is larger than that for pi/ 2 states. Moreover, using Fig. [7|and 
Eqs. JSJ) we conclude that at large Z the finite nuclear size effect on Lamb shift of P3/2 states is only by a few times 
smaller than that for s states. These conclusions are not consistent with numerical data at large Z . Certainly, our 
calculation for p 3 / 2 states is valid only up to the order a(Za). The p 3 / 2 correction does not contain the ultraviolet 
logarithm and therefore the convergence of Za-expansion can be different from that for s\/ 2 and p\/ 2 states. A poor 
convergence of the expansion could be the reason for disagreement with previous numerical calculatins. However a 
new exact in Za numerical calculation for p 3 / 2 states would be interesting. 

We are grateful to G. Soff for helpful comments. The work was supported by RFBR grants 01-02-16926 and 
03-02-16510. 
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